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Polycrystalline samples of the recently discovered MgCNis superconductor were investigated by 
transport, ac susceptibility, dc magnetization and specific heat measurements in magnetic fields up 
to 16T. Consistent results were obtained for the temperature dependence of the upper critical field 
Hc2 [T) from resistance, ac susceptibility and specific heat measurements. A WHH like temperature 
dependence of Hc2 (T) and the quadratic relationship Hc2 (0) ~ point to an effective predominant 
single band behavior near the quasi clean limit. Evidence for strong electron-phonon and electron- 
paramagnon coupling was found analyzing the specific heat data. The corresponding s- and p-wave 
scenarios are briefly discussed using calculated densities of states of different Fermi surface sheets. 
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I. INTRODUCTION 

The recent discovery of superconductivity iii the inter- 
metaUic antiperovskite compound MgCNialll with a su- 
perconducting transition temperature of Tc — 8 K is 
rather surprising considering its high Ni content. There- 
fore, it is not unexpected that this compound is near a 
ferromagnetic instabiHtyjwhich might be reached by hole 
doping at the Mg sites .□ The possibiHty of unconven- 
tional superconductivity due to the proximity of these 
two types of collective order has attracted great interest 
in the electronic structure and the physies,ftf the pairing 
mechanism. Band structure calculationsucHa for MgCNia 
reveal a domination of the electronic states at the Fermi 
surface by the 3d orbitals of Ni. 

MgCNis can be considered as the 3-dimensional ana- 
logue of the layered transition metal borocarbides which 
have superconducting transition temperatures up to ~ 
23 K. In spite of the much lower Tc of MgCNia, its 
upper critical fields rSinSn^t low temperatures is, with 
Hc2(0) = 10 ... 15 TBOQ'I comparable with that of the 
borocarbides or even higher. This is connected with the 
completely different temperature dependence of ffc2 for 
these compounds. The Hc2{T) dependence of MgCNia is 
similar to that of dirty-limit intermetallic superconduc- 
tors with a steep slope of H^2{T) at T^. 

Through analysis of specific heat data, MgCNia was char- 
acterized in the framework of a comjentional, phonon^ 
mediated pairing both as a moderatetlB and as a strongLl 
coupHng superconductor. Strong coupHng is also sug- 
gested by the large, energy gap determined from tun- 
neling experiments.!] The question of the pairing sym- 
metry is controvepfiially discussed in the literature. ^"^C 
NMR experimental support s-wave pairing in MgCNis, 
whereas tunneling spectra indicate an unconventional 



pairing state.Ll 

In the present investigation, upper critical field and spe- 
cific heat data of MgCNia were analyzed with the aid of 
theoretical results for Fermi velocities and partial densi- 
ties of states in order to find out a consistent physical 
picture explaining the experimental results. 



II. EXPERIMENTAL 

PolycrystaUine samples of MgCNis have been prepared 
by solid state reaction. In order to obtain samples with 
high Tc; one has to use an excess of carbon as stated 
m Ref. 0. To cover the high volatility of Mg during 
the [Sintering of the samples an excess of Mg is needed, 
too.EJ In this study, a sample with the nominal formula 
Mg1.2C1.6Ni3 has been investigated which is denoted as 
MgCi.gNis. To prepare the sample, a mixture of Mg, 
C and Ni powders was pressed into a pellet. The pellet 
was wrapped in a Ta foil and sealed in a quartz ampoule 
containing an Ar atmosphere at 180 mbar. The sample 
was sintered for half an hour at 600° C followed by one 
hour at 900° C. After a cooling process the sample was 
regrounded. This procedure was repeated two times in 
order to lower a possible impurity phase content. The 
obtained sample was investigated by x-ray diffractome- 
try to estimate its quality. The diffractometer pattern 
(Fig. 0) showed small impurity concentrations mainly re- 
sulting from MgO and unreacted carbon crystaUized in 
form of graphite. The lattice constant of the prepared 
sample was determined to be a = 0.38107(lljjim using 
the Rietveld computer program FULLPROFJij Accord- 
ing to Ref. ^ this indicates that the nearly single phase 
sample corresponds to the superconducting modification 
of MgCxNis. 
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Figure 1: Rietveld refinement for MgCi.eNia. The crosses 
correspond to the experimental data. The solid line shows 
the calculated pattern. Vertical bars give the Bragg positions 
for the main phase MgCNia, graphite and MgO (from top to 
bottom). The line at the bottom gives the difference between 
the experimental and calculated data. 



The superconducting transition of the sample was inves- 
tigated by measurements of the electrical resistance, the 
ac susceptibility and the specific heat. For the electri- 
cal resistance measurement a piece cut from the initially 
prepared pellet with 5mm in length and a cross section 
of approximately Imm^ was measured in magnetic fields 
up to 16 T using the standard four probe method with 
current densities between 0.2 and 1 A/cm^. The ac sus- 
ceptibility and the specific heat measurements were per- 
formed on other pieces from the same pellet in magnetic 
fields up to 9 T. 



III. RESULTS 
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Figure 2: Temperature dependence of resistance of MgCi.eNis 
at zero applied magnetic field. The inset shows the supercon- 
ducting transition region. 
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Figure 3: Field dependence of resistance of MgCi.eNis mea- 
sured at several temperatures. 



A. Superconducting transition and upper critical 
field 

In Fig. |[ the temperature dependence of the electrical 
resistance of the investigated sample is shown. A super- 
conducting transition with an onset (midpoint) value of 
Tc = 7.0 K (6.9 K) is observed (see inset of Fig. |^) which 
is consistent with the onset of the superconducting tran- 
sition at Tc = 7.0 K determined from ac susceptibility. 
It should be noted that the sample shown in Fig. || has a 
resistivity of paooK = 2.1 milcm at 300 K which is much 
too large in order to be intrinsic. On the other hand, its 
residual resistance ratio RRR = i?(300K)/i?(8K) = 1.85 
and the shape of the R{T) curve are typical for MgCNia 
samples.EJ A possible explanation for the high resistivity 
of the investigated sample which was not subjected to 
high pressure sintering is a relatively large resistance of 
the grain boundaries. 



The field dependence of the electrical resistance of the 
same sample is shown in Fig. || for several temperatures 
between 6.0 and 1.9K. A sharp transition is observed 
which remains almost unchanged down to low temper- 
atures. In Fig. 1^, the fields -ffio, H^o and Hgo defined 
at 10%, 50% and 90% of the normal-state resistance are 
plotted as function of temperature. Identical results have 
been found from resistance-vs.-temperature transition 
curves measured at different magnetic fields. Addition- 
ally, Fig. ^ shows upper critical field data determined 
from ac susceptibility measurements. The onset of su- 
perconductivity was used to define Hc2 from ac suscep- 
tibility. 

It is clearly seen that for the investigated sample H^2^ 
{Hc2 obtained from susceptibility) agrees approximately 
with HiQ. A similar behavior was already observed for 
MgB2, whereas in the case of rare-earth nickel borocar- 
bides the onset of superconductivity determined from ac 
susceptibility was typically found to agree well with the 
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Figure 4: Temperature dependence of the upper critical field 
for MgCi.eNis. The circles show the midpoint value -ffso of the 
resistivity in the normal state. The two lines labelled Hio and 
Hso denote 10% and 90% of the normal state resistivity. The 
triangles represent the upper critical field determined from 
the onset of ac susceptibility. 



Figure 5: Upper critical field at zero temperature vs. super- 
conducting transition temperature for the MgCi.eNis sample 
of this work ( ■ ) and of several MgCNis samples reported in 
Ref. 3 ( O ), Ref. 4 ( • ), Ref. 5 ( ▼ ) and Ref. 6(a). The 
experimental data can be described by a quadratic law (solid 
line). 



midpoint value (i?5o) of the normal state resistivity. The 
width AH = Hqq—Hiq of the superconducting transition 
curves in Fig. |^ (and Fig. ^ remains, with AH ~ 0.6 T, 
almost unchanged down to low temperatures. A poly- 
crystalline sample of a strongly anisotropic superconduc- 
tor shows a gradual broadening of the superconducting 
transition with decreasing temperature as was observed, 
for example, for MgB2.Ej Therefore, the nearly constant 
transition width AH observed for the investigated sam- 
ple can be considered as an indication of a rather small 
anisotropy of Hc2{T) in MgCNia. 

The extrapolation of HgoiT) to T = yields an upper 
critical field of i?c2(0) — 11-3 T. The observed tempera- 
ture dependence of the upper critical field is typicaLfMj 
Hc2{T) data reported for MgCNiai-ajnd was describedOd'H 
within the standard WHH modeta by conventional su- 
perconductivity in the dirty limit. The WHH model pre- 
dicts a relation 77c2(0) oc {—dHc2/ dT) -Tg. Available 
data for MgCNia, including those presented in this paper 
show a strong variation of i?c2(0) with Tc (see Fig. 5), 
whereas [dH^2/dT)^^ « -(2.65 ± 0.2) T/K remains al- 
most unchanged. Considering the data in Fig. 5, the 
linear dependence of i?c2(0) on Tc predicted for constant 
{dH^2/ dT)rp by the WHH model can be ruled out. 

The solid line in Fig. || corresponds to a quadratic law 
Hc2{^) ~ which is the benchmark for the clean limit. 
IndjSfid in the isotropic single band s-wave clean limit one 
hastil 

i?c2(0) [Tesla] = 0-0237 ^^g^ (1) 

Compared with WHH the effect of strong coupling (mea- 
sured by the dimensionless electron-phonon coupling con- 



stant A) is an enhancement of i?c2 primarily through 
the renormalization of the bare Fermi velocity up — > 
wf/(1 + A) and a weak enhancement factor coming from 
the energy dependence of the gap function. The strong 
coupHng (A « 3.4 near the Brillouin zone boundary and 
zero at the zone center) to the rotational mode near 13 
meV proposed in Ref. ^ gives an upper limit for a strong 
coupling required to reproduce upper critical fields as 
high as 12 ... 15 T. Using the averaged Fermi velocities 
of 1.45 • 10^ m/s for the rounded cube-like Fermi surface 
shown in Fig. O we arrive at A w 2.5 ... 3. 



B. Specific heat 

Specific heat measurements were performed in order to 
get additional information about Hc2{T), the super- 
conducting pairing symmetry and the strength of the 
electron-phonon coupling from thermodynamic data. In 
Fig. ^ specific heat data, Cp/T vs. T^, are shown for 
applied magnetic fields up to 8T. The upper critical 
fields, Hc2(r), determined from the specific heat data, 
are shown in Fig. |^. It is clearly seen that the Hc2{T) 
data obtained from the specific heat are located in the 
small field range between the Hgo(T) and Hiq{T) curves. 

In order to describe the experimental data of the normal 
state specific heat in zero field between and ZQK (see 
Fig. H) the expression 

Cn{T) — 7nT + Clattice + CEinstein + J^CSchottky (2) 

was used. The linear-in-T term is due to the electronic 
contribution with 7n as the Sommerfeld parameter. The 
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Figure 6: Specific heat data Cp/T vs. of MgCi.eNis mea- 
sured at various magnetic fields up to 8T. The solid line is a 
fit of Eq.(^) to the data for H — above Tc. Its intersection 
with the ordinate gives the Sommerfeld parameter 7n. 
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Figure 7: Comparison of upper critical field data determined 
from specific heat ( • ) and resistance measurements for 
MgCi.eNia. Hw and Hgo were determined at 10% and 90% 
of the normal state resistivity, respectively. An entropy con- 
serving construction was used to determine the upper critical 
field from the specific heat data of Fig. B. 



second term, ciattice = PT^ + <5T^ represents the phonon 
contribution. This extension of the usual Debye approx- 
imation Ciattice oc includes deviations from the linear 
dispersion of the acoustic modes and is required in or- 
der to describe the phonon contribution to the specific 
heat in an extended temperature range. The same pro- 
cedure ha&,already been used to describe the specific heat 
of MgB2]i3 An additional single Einstein mode, CEinsteim 
was used to take into account the contribution of the 
previously mentioned lowest energy optical mode near 
13meV.a The last part, ?^Cschottky, represents the contri- 
bution of a conventional two-level Schottky model which 



was used to describe the slight upturn of the specific heat 
just above the superconducting transition temperature 
(see Fig. ||). This feature possibly results from contribu- 
tions of paramagnetic impurities like unreacted Ni (the 
prefactor n gives the concentration of these impurities) . 

The fitting procedure returned a Sommerfeld parameter 
7n = 35mJ/niplK^ which is close tp the reported values 
of 7n = 27.m and 33mJ/molK2.y The Debye temper- 
ature was found to be 6d = 289K and the energy of 
the optical phonon mode resulting from the fit is 9.5meV 
which slightly deviates from the theoretically proposed 
value (13meV)fl The paramagnetic impurity concentra- 
tion derived from the fit is n « 4.5% and the energy gap 
of the Schottky anomaly is 5.3meV. As shown in Fig.Js], 
the experimental data are very well described by Eq. (Q) 
in the investigated temperature range < T < 30K. 
The jump Ac of the specific heat at (see inset of Fig. ||) 
is given by the difference between the experimental data, 
Cp and the normal specific heat contribution, Cn- 

The transition temperature = 6.83K obtained from 
the inset of Fig. ^ agrees approximately with the tran- 
sition temperatures = 7. OK and = 6.9K derived 
from ac susceptibility and from resistance data, respec- 
tively. In the low temperature region the experimental 
data for Ac/T versus T can be described by the BCS-like 
expression 



Ac = 7.957nTc exp 



A(0) \ 

kBTj 



(3) 



using 2A(0)/fcBl'c = 2.96 instead of 2ABcs(0)/fcBTc = 
3.52 predicted by the BCS model. Above T = 4K the 
fit of Eq. (^ starts to deviate from the experimental 
curve resulting in a much lower value for the jump at 
T — Tc than observed experimentally. Notice that the 
experimental value of the jump, Ac/7nT'c = 1-6, nearly 
corresponds to the BCS value Ac/7nT'c = 1-43 indicating 
weak electron-phonon coupling. This contrasts to the 
strong coupling derived from Eq. (|l|) and Hc2{0) data. 
A natural explanation for this discrepancy is the two- 
band character of MgCNia which will be discussed in the 
next section. 

To examine the temperature dependence of the electronic 
specific heat 



Cei(r) = Ac + 7nT 



(4) 



at iJ = in detail, Cei(T) / 'fj^Tc is plotted logarithmi- 
cally vs. Tc/T in Fig. It is clearly seen, that the 
experimental data at low temperatures {T^/T > 2) fol- 
low the modified BCS expression (Eq. d), solid line). 
We found that this exponential law is not affected by the 
low temperature branch of the Schottky term in Eq. (||) 
which has an exponential temperature dependence, too. 
The exponential temperature behavior of the electronic 
specific heat at low temperatures found for MgCNia is 
a strong indication for s-wave superconductivity in this 
compound. 
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Figure 8: Specific heat data Cp/T vs. of MgCi.eNis at zero magnetic field. According to Eq.(|2|), the dashed line represents the 
lattice contribution, the dotted line shows the contribution from a rotational phonon mode (see text) and the dash-dotted line 
gives a Schottky contribution. Inset: Contribution of the superconducting electrons to the specific heat Ac/T = (cp — Cn) /T. 
The solid line in the inset is a fit of the BCS expression to the data with parameters according to Eq. (^). The conservation 
of entropy was confirmed by integrating of Ac/T in the temperature range < T < Tc (according to the data above 2K and 
the solid line below 2K). 




Figure 9: Normalized electronic specific heat contribution of 
MgCi.eNia vs. Tc/T. The comparison of the data with the 
BCS expression (solid line) corresponding to Eq. (jsj) clearly 
shows the exponential temperature dependence of the elec- 
tronic specific heat at low temperatures. 



In the superconducting state, a linear-in-T electronic 
specific heat contribution ^{H)T arises from the nor- 
mal conducting cores of the flux lines for applied mag- 
netic flelds H > H^i. This contribiition can be ex- 
pressed as "f{H)T — Cp (T, H) — Cp (T, 0) ,t3 where Cp (T, 0) 



is the speciflc heat in the Meissner state. Speciflc heat 
data for MgCi.eNis at T = 2K (see inset of Fig. 10) 
were analyzed in order to derive the field dependence 
of l{H). In Fig. the obtained 7(i?)/7N is plot- 
ted against H/Hc2{0) using the Sommerfeld parameter 
7n = 35mJ/molK2 and i?c2(0) = 11.3T. Not shown in 
this plot are high-field data which are infiuenced by an 
additional contribution to the specific heat in the normal 
state arising in magnetic fields. This contribution which 
is clearly seen in Fig. |^ as deviation of the Cp/T data 
from the solid fine causes a shift of the high-field Cp/T 
data in the superconducting state to higher values. The 
origin of this effect is not yet understood. Its infiuence 
on the field dependence of Cp/T at T = 2K is illustrated 
in the inset of Fig. 10. At low magnetic fields, a negative 
curvature is observed which starts to change its sign at 
fields above 4T. 

The low-field data of Cp/T shown in Fig. |l^ can be de- 
scribed by the expression 7/7N = (i?/-ffc2(0)) ' which 
differs from the Hnear j{H) law expected for isotropic 
s-wave superconductors in the dirty limit. A non-linear 
field dependence close to j{H) cx has been reported 
for some unconventional superconductors with gap nodes 
in the quasiparticle spectrum of the vortex state as in 
YBa2Cu307jlj and in the heavy fermion superconduc- 
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Figure 10: Field dependence of the specific heat contribution 
'y{H) of the vortex core electrons in the mixed state normal- 
ized by the Sommerfeld parameter 7n and -ffc2(0). The black 
line is a fit of 7(-ff)/7N = {H/ Hc2{0))"'^ to the experimental 
data. Inset: Field dependence of the specific heat Cp/T at 
T = 2K. 



n s-wave superconduc- 
and the borocarbides 



tor UPtsjEJ but also-in some|d|6 
tors such as CeRu2,B NbSjeepI 
RNi2B2C (R = Y, Lu)Eill23 DelocaHzed quasiparticle 
states around the vortex cores, similar to these in d-wave 
superconductors, seem to be responsible fos- 
linear 7(i?) dependence in the borocarbidesB^ 



IV. DISCUSSION 



In principle, the upper critical field data found for 
MgCNia could be understood both in the s- and p-wave 
scenario. The high magnitude of the upper critical fields 
might be achieved in s-wave superconductors, but also in 
a clean limit weak coupling p-wave case employing the 
"slow" holes on the "four-leafed clover"-like Fermi surface 
sheets-jwith 0.5 • lO^m/s (see Fig. [u]). Following Maki 
et al.E^ an additional numerical factor 1.3 should be in- 
troduced in this case in Eq. (0). The different temper- 
ature dependence of the electronic specific heat at low 
temperatures in s- and p-wave superconductors allows 
the discrimination between predominant s- and p-wave 
scenarios. As mentioned above, the exponential temper- 
ature behavior of the electronic specific heat found for 
MgCNis is a strong indication for s-wave superconduc- 
tivity in this compound. 

It is convenient to rewrite Eq. (Q) using experimentally 
accessible quantities as the plasma energy Wpi, the vol- 
ume of the unit cell V, and the Sommerfeld constant jn- 
Then, Eq. (||) provides a criterion for a superconduc- 
tor which can be described in the clean limit within the 



isotropic single band model: 

3.6c^2^ [eV2] i/e2(0) [T2] V [(10-iOm)' 



Q = 



7n [mJ/molK2] [K] (1 + A) 



1.4 



(5) 



If Q differs significantly from 1, a more complex model 
should be considered. 

Some values for Q are summarized in Tab. |. 
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Table I: Values of Q and further parameters (see text) for se- 
lected superconductors. Q is a measure for the applicability 
of the isotropic single band model. A deviation from Q ^ 1 
indicates the need for a more complex model. For MgB2 and 
MgCNis, two limiting values for A are consiered for illustra- 
tion. 

The crystal structure of MgCNis can be seen as a 
three dimensional analogue of the layered borocarbides. 
Thus it is instructive to compare the results with the 
electronic specific heat dependence of the borocarbides 
and the boronitrides which were found to be isostruc- 
tural to the borocarbides. The electronic specific heat 
of La3Ni2B|gJ^l3_5 exhibits an exponential temperature 
dependenceEa as does MgCNis, while that of YNi2B2C 
or LuNi2B2C follows a power law of the type Cei = 

37n (t~) ■^ith an exponent a ~ 3.i3S 

The electron-phonon coupling constant A in MgCNia av- 
eraged over all Fermi surface sheets can be estimated 
from the relation 



7n 



-4(l + A)iV(^F) =7o(l + A) 



(6) 



using the experimental value of the Sommerfeld constant 
7n = 35mJ/molK2 and the density of states (DOS) 
at the Fermi level N{Ey) calculated within the local 
density approximation. From the calculated density of 
states and the corresponding bare specific heat coeffi- 
cient, 70 = llmJ/molK^ one obtains a large electron- 
phonon coupling constant A = 2.2. 

According to band structure calculations using the 
full-potential nonorthogonal local-orbital minimum-basis 
schemecj within the local density approximation, the to- 
tal DOS can be decomposed into a roughly 15% and a 
85% contribution stemming from the fast and the slow 
sheets of the Fermi surface, respectively. Then, the total 
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Figure 11: The two Fermi surface sheets of MgCNis and the 
corresponding band resolved density of states near the Fermi 
level. "Band 1" corresponds to the Fermi surface sheet in the 
topright panel, "band 2" to the Fermi surface in the topleft 
panel. 



coupling constant averaged over all Fermi surface sheets 
reads 



At, 



Ni(0) N2(0) 

Ai h Ai- 



N(0) 



N(0) 



(7) 



and with the aid of Atot = 2.2 and Ai w 2.5 ... 3 esti- 
mated from the upper critical field i?c2(0) ~ 12 . . . 15T, 
respectively (Eq. (gj)) one arrives at A2 ~ 2.0 . . . 2.1. The 
question then arises about the origin of these strong mass 
enhancements in both bands. Obviously A2 cannot be of 
electron-phonon nature alone, since otherwise one should 
ask, why are Tc, the above mentioned gap, and Hc2{0) so 
low? So we are forced to assume that A2 should be decom- 
posed into an electron-phonon contribution and a pair- 
breaking electron-paramagnon As/ (electron-spin fluctu- 
ation) one. We adopt for the electron-phonon part a typ- 
ical transition metal value, say X2ph 1. From Fig. ^ a 
band width of about 0.6 eV can be estimated for band 2. 
Taking the average of the electron and the hole Fermi en- 
ergies as a representative effective Fermi energy of band 
2 hw2F ~ 0.3 eV can be estimated. Then, according to 
the Berk-Schrieffer theory the paramagnon spectral den- 
sity should exhibit a maximum near hugf ~ lioj2F/S ~ 
80 meV, where S" « 4 is the Stoner factor. Since this 
frequency exceeds considerably the typical phonon fre- 
quencies, the effect of the paramagnon pair breaking is 



somewhat reduced due to the pseudopotential effectEa, 
well-known for the large bare Coulomb repulsion /x where 
it produces the smaller pseudopotential ^* -value due to 
the large logarithmic factor of about 5 entering its de- 
nominator. In the present case, this factor is reduced to 
1.4, only. Then the "bare" paramagnon coupHng constant 
Xsf 1 is reduced to about A*^ ~ 0.4. Together with 
standard Coulomb pseudopotential, e.g. 0.13, one never- 
theless arrives at a sizable suppressed pairing in band 2 
which is responsible for the observed small gap seen in 
the low T specific heat data. The missing od only very 
weakly pronounced curvature in Hc2{T) near Tc points to 
weak interband coupling, a posteriori justifying our one 
band estimate for i?c2(0) given above. A more quanti- 
tative study must await more detailed knowledge on the 
phonon and paramagnon spectra. 



CONCLUSION 



The electronic specific heat shows an exponential temper- 
ature dependence at low temperatures which is a strong 
indication for s-wave pairing symmetry in this compound. 
Strong electron-phonon coupling has been derived from 
specific heat data and the calculated densitiy of states av- 
eraged over all Fermi surface sheets. Comparing MgCNia 
with other multi-band superconductors such as MgB2 
or transition metal borocarbides, one concludes that the 
disjoint Fermi surface sheets differ from each other not 
only in the strength of the electron-phonon interaction 
or the Fermi velocities but to the best of our knowledge 
also in strength of the depairing interaction. Approach- 
ing the remarkable peak in the DOS slightly below the 
Fermi energy, an interesting enhancement of the param- 
agnon contribution followed by a gapless Fermi surface 
sheet and a possible p-wave superconductivity mediated 
by ferromagnetic spin fluctuations might be expected. 
In any case, the competing interplay of strong electron- 
phonon and electron-paramagnon interactions together 
with a strongly energy dependent DOS on disjoint, al- 
most decoupled Fermi surface sheets is a great challenge 
for future theoretical and experimental work. 
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